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Abstract 



We present observations of the full January 2010 transit of HD80606b from 
the Canadian microsatellite, Microvariability and Oscillations of Stars (MOST). 
By employing a space-based telescope, we monitor the entire transit thus limit- 
ing systematic errors that result from ground observations. We determine mea- 
surements for the planetary radius (R p = 0.987 ± 0.061Rj up ) and inclination 
(i = 89.283° ± 0.024) by constraining our fits with the observed parameters of 
different groups. Our measured mid-transit time of 2455210. 6449±0. 0034 HJD is 
consistant with the 2010 Spitzer results and is 20 minutes earlier than predicted 
by groups who observed the June 2009 transit. 

Subject headings: transits — planets and satellites: individual (HD80606b) - 
techniques: photometric 
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Introduction 



With its unusually high eccentricity of e = 0.93 and long orbital period of 111 days, 
HD80606b challe nges our understanding of exoplanets. This planet was discovered by 



NaefetaL 



( 1200 ll ) via radial velocity techniques in which they determined a minimum 
planetary mass of 3.9 Mj up . The occultation of HD80606b was observed with the Spitzer 
Space Te lescope eight years la ter indicating that the inclination of the planet is greater 



than 89° (lLaughlin et al. 



2009). The existence of an occultation indicated that there was a 



15% chance that this planet would also have a visible primary transit. 

Despite the low probability, several independent groups detected the Febr u- 



ary 2009 primary transit from ground-based observat 


ories ( 


Fossev et al. 


2009 




Garcia-Melendo & McCulloueh 2009: Hidas et al. 


2010; 


Pon 


t et al. 


2009 


; I 


Vloutou et al. 


2009). The full transit was later observed by 


Winn et al. 


f 


2009 


) who used multiple 



ground-based observatories to sew together complete coverage. 



predicted by 



Winn et al.l (120091) also confi rmed a spin-orbit misalignment of 53° that was first 



Moutou et al. 



(120091 ). It is possible that the high eccentricity and tilted orbit 



of HD80606b are the result of the Kozai m echanism induced by the system's binary stellar 



companion HD80607 (IWu fc Murray 



2003|). 



Since the 12 hour dur ation of HD806 6b's t ransit presents a difficulty for accurate 



ground-based observations, 



Hebrard et al 



(120 101 ) employed the combination of Spitzer 



photom etry and SOPHIE sp ectroscopy to continuously observe the full January 2010 



transit. 



Hebrard et al. 



(120101 ) discovered t heir mid-transit tim e to be appr oximately 24 



minut es earlier than the value predicted by 



Winn et al. 



(120091 ). In addition, 



Hebrard et al 



Pont et al. 


(2009 


) and 


Winn et al. 


(2009) 
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Shporer et al.l (120101 ) launched a 



multiple-site g r ound- based observation of the same 



transit in an attempt to confirm the 



Hebrard et al 



(120101 ) timing discrep ancy. Their 



midtransit time was 12 minutes earlier than predicted and within 1.3a of 



Hebrard et al. 



( I2010[ ). 



We observe the January 2010 transit of HD80606b with the Microvariability and 



Oscillations of Stars (MOST) earth-orbiti ng optical te l escop e wit h the goals of: 



determining a mid-transit time using both 



Winn et al. 



(120091 ) and 



Hebrard et al 



01 



(120101 ) 



orbital parameters and (2) derive an independent set of results for the planetary radius, 
inclination, and mid-transit time. Our observations by the MOST telescope are described 
in §2, our light curve modeling is detailed in §3, and we present and discuss our results to 
both goals mentioned above in §4. We conclude this paper in §5. 



2. Observations 



The Microvariability and Oscillations of Stars (MOST) is a microsatellite with a 15 



cm aperture photometer designed for aster oseismology, 



3ut is well-suited to observe short 



period exoplanets around solar type stars (iRowe et al 



20061). The satellite employs a 



broadband visible filter centered at 525 nm. MOST is in a low polar, sun-synchronous orbit 
that permits for a continuous viewing zone (CVZ) of stars between the declination of -19° 
and +36°. HD80606b, however, has a declination of +50° placing it outside of the CVZ. 
Therefore, we experience short blackouts in our data as Earth passes in between MOST and 
HD80606. We are unable to observe the January 8th, 2010 secondary eclipse of HD80606b 
due to an inconvenient instrumentation issue that led to a 12 hour gap in our data. 

We observe HD80606b starting from 24548677.97302083 HJD and continuing to 
2454877.85480324 HJD (January 7th to January 19th, 2010) with an exposure time of 80.40 
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second integrations. Since HD80606 is a magnitude 9.06 star, we use the direct imaging 
mode of the CCD and employ simple aperture photometry with a field of view of 20 x 20 
pixels. We capture images of both HD80606 and HD80607, but we are able to constrain 
HD80606 within a 3 pixel radius compared to the binary stars' separation of 10 pixels. 

Stray light for MOST varies from 100 ADU to 1500 ADU as it passes over the polar 
caps and day/night boundary of Earth. We remove data that has a background level above 
800 ADU (about 100,000 points). At this lower level of stray light, our photometric error is 
equal to the expected Poisson noise li mit. For a de t ailed procedure of standard observation 



and data reduction of MOST data, see 



RoweetaL 



(I2006h 



Figure [T] shows our full photometric time series, with points binned every 15 minutes 
for clarity. By binning the data the gaps are not as apparent; however, we fit the full 
unbinned transit (Section 3). Due to the binning effect, the transit in Figure [1] appears 
asymmetric. A close up of the unbinned transit in Figure [2] demonstrates that within the 
error of our measurement, the transit is symmetric and we are not forced to make necessary 
modifications for an asymmetric effect. The large gap occurring at 24552005.5 is due to the 
instrumentational error discussed above. 



3. Lightcurve Modeling 



We use the fitting process described in 



Barnes fc Fortneyl (120031 ) . With high eccentricity, 



the transit duration is correlated to the inclination and argument of periapsis of the planet's 
orbit. The transit ingress and egress times also experience a several minute asymmetry, 
with the ingress duration shorter than the egress in the case of HD80606b. We account for 
these y a riatio ns by explicitly incorporating orbital eccentricity in to our code, as detailed in 
Barnesl fj2007n . 
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In this routine, we do not account for the impact parameter b as this parameter is 
directly related to the inclination in an eccentric orbit by 



a 



1 



COS[l) 



1 + e * sin{u)) 

By fitting for the inclination alone, we accelerate the fitting algorithm. 



(1) 



We fit the out-of-transit flux dynamically in our model instead of assuming a value. 
Out-of-transit points that fall more than 3a outside of this photometric baseline are removed. 
Our best fit transit model is determined by minimizingx 2 through a Levenberg-Marquardt 
fitting algorithm. 

As our data is not precise enough to fit for the limb darkening parameters combined 
with our other fitted parameters, we hold both coefficients constant. Our limb darkening 
parameters are established by using th e defined MOST quadratic parameters of HD80606 



with a T e ff = 5QA5K and log g = 4.50 (iNaef et al 



200JJ). These co efficients are manip ulated 



to fit our model by using the limb darkening method highlighted in 



Brown et al. 



(120011 ). Our 



first coefficient, designated as C\ = 0.742, represents the magnitude of d arkening while our 



secon d coefficient, C2 = 0.458, accounts for the curvature of the transit (IBarnes Sz Fortney 



2003 1. 



We do not account for any other systematic effects or correl ated noise. Therefore, we 



assume all err ors to be random with a Gaussian distribution. See lBarnes &: Fortneyl ( 120031 ): 
Barnes! ( 120071 ) for additional details regarding our model. 



We execute t wo different m o del fi ts: the first model is used to verify the mid-transit 
time reported by 



Hebrard et al 



(120101 ) an d a second to determ ine our own independent 



parameters. In our effort to reproduce the 



Hebrard et al. 



(120101 ) mid-transit time, we hold 



all their parameters consta nt and fit so l ely fo r the baseline flux and mid-transit time. We 



repeat this procedure using 



Winn et al 



( 120091 ) parameters for comparison to their predicted 
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mid-tr ansit time. For clarity , we will den ote our timing resu lts with either a superscript 



TH for 



Hebrard et al. 



(120101 ) and TW for 



Winn et al 



(120091 ) when discussing out results in 



Section O 

In establishing our independent best-fit parameters, we determine that our data are 
only adequate to fit for four parameters: planetary radius (R p ), mid-transit time (T m id), 
incli nation (i), and the ba selin e flux. We constra in our two fits with the parameters 



from 



Hebrard et al. 



(bold ) and 



Winn et al 



(120091 ) for stellar radius, stellar mass, period, 
eccentricity, argument of periapsis, and semi-major axis. We analyze these results in Section 
E2J 



4. Discussion 

4.1. Mid- Transit Time Fit Only 

In this section, we discuss our mid-transit time results of HD80606b with a goal 
of verifying the early mid-transit time discovered by 



Hebrard et al. 



(2010). In order to 



measure a mid-transit time appropriate for comparison, we hold iHebrard et al. 



( 20101 ) 



orbital, stellar, and planetary parameters constant and fit only for the mid-transit time and 
baseline flux. 

Our resulting mid-transit t ime (d enoted with the su perscript TH) is 20 m inutes earlier 



than predicted by 



Winn et al. 



(120091 ) . For comparison, 



Hebrard et al. 



(I2010f ) 



mid-tr ansit time 24 minutes (17a difference) earlier than the value predicted by 



bund a 



Winn et al. 



(120091 ) . 



Orbital paramet ers have a direct effec t on mid-transit times, especially in a highly 

I2OIO ) determined a separate set of orbital parameters 



compared to 



t. Since 


Hel 


Drard et al. 


Winn et al. 


( 


2009 


). we 



(120091 ) . we also fit the MOST lightcurve using the latter's 
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parameters for the mid-transit time. Even with the different orbital parameters, we 
determine a mid-transit time (denoted with the superscript TW) 20 minutes early as well. 

We conclude that the difference in mid-transit times between the two groups is 
not dependent on their varying orbital parameters. Furthermore, we determine that 
our results have less than a la discrepancy from the mid-transit time determined by 



Hebrard et al. 



(120101 ). A visual comparison of each mid-transit time with uncertainties is 
displayed in Figure E l We also high li ght th e mid-transit timing results from an independent 



ground-based group, IShporer et al. 



(120101 ) . in Figure [3] who observed the 2010 January 



transit. Our overall timing results are also listed in Table [TJ 



(120101 ) result to a less than la 



Our mid-transit time agrees with the iHebrard et al. 
difference. We also find that our m id transit time has l.Qa difference with the result 

(2010) and a greater than 3cr difference from the predicted 



determined by 
transit time of 



Shp orer et a 



Winn et al 



)9) 



(120091 ) . I n fact, the only resu 



(120091 ) predicted time (within la) is 



Shporer et al. 



t that agrees with the 



Winn et al 



(120101 ). They note, however, that due to 



the difficulty o 



Shporer et al. 



grou nd based observations, the uncertainty in their measure ment is large . 



(120101 ) still find a mid-transit time of 12 minutes earlier than 



Winn et al. 



(120091 ) predicted, a large discrepancy regardless of the uncertainty. 



One possi ble explanat i on for the large difference between the mid-transit times of 



this work and IWinn et al 



(120091 ) is the possibility of a perturbing bo dy in the system 



creating a transit timing variation (TTV). However, 



another body pertur 



NaefetaL 



(120011 ) show that if there is 



Ding this system, it is not detectable via radial velocity measurements. 
Wu fc Murray! (J2003|) also argue that while another planet may have formed around 
HD80606, no Jupiter-mass planet could lie in a stable orbit between 0.05 to 100 AU. They 
further show that no Earth-mass planet could exist between 1 and 20 AU. As TTVs are 
sensitive even to small mass objects, there is the possibility of the existence of a moon 
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around HD80606b. Contrary to this, iBarnes &: O'Brienl (120021 ) show that no moon massive 
enough to create a transit timing variation could exist in a stable orbit around the highly 
eccentric planet. 

With the combination of 5 different mid-transit timing results, the possible mid-transit 
time of HD80606b spans over 24 minutes. The large discrepancy between each group's 
result demonstrates the mid-transit timing issue is unlikely due to an astrophysi cal 



Hebrard et al 



(|2010h (both 



event. Furthermore, it is of interest that both this work and 
space-based observations) agree within less than la difference of an early mid-transit 
time, while the ground-based results determine a mid-transit time several minutes later. 
We therefore propose the hypothesis that the difference in mid-transit time derives from 
systematic errors that depend on the method used to observe HD80606b. 





Mid- Transit Time 




MOST (This Work) TH 


2455210.6449 ± 0.0034# JD 


1.11735 


Spitzer (Hebrard et al. 2010) 


2455210.6420 ± 0.0010#JX> 




Ground-Based (Shporer et al. 2010) 


2455210.6502 ± 0.0064# JD 




MOST (This Work) TW 


2455210.6447 ± 0.0032# JD 


1.11436 


Ground-Based (Winn et al. 2009*) 


2455210.6590 ± 0.0050# JD 





*Predicted 

Table 1: Our mid-transit times for the January 2010 transit. "MOST Ti? " designates mid- 
transit time using parameters determined by Hebrard et al. 2010. "MOST TW/ " represents 
the mid-transit time determined with Winn et al. 2009 parameters. Xr 2 shows the reduced 
X 2 of our best fit. 
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4.2. Complete Transit Fit Parameters 



We determine two independent sets of parameters for HD80606b, which we will discuss 
and compare in this section. Figure [2] shows MOST's HD80606b transit lightcurve with the 
data unbinned and each point representing 80.40 seconds of observation. The large gaps 
between each cluste r of data is the result of Earth passing between MOST and HD80606 



(I Walker et al. 



20031 ) . Unfortunately, two of these blackouts occur during large parts of the 
transit ingress and egress. In order to minimize the effect of these blackouts on our fit, we 
hold our limb darkening coefficients constant during the entire fitting routine. 

We further minimize the effects of scatter and blackouts in our data by electing to fit 
for only four parameters. These four parameters are: planetary radius (R p ), inclination (i), 
mid-transit time (T m id), and the out-of-transit baseli ne flux. In this sectio n, we will denot e 



our re sults with either a superscript FH or FW, for 



Hebrard et al. 



(I2010h or 



Winn et al. 



( 120091 ) respectively, based on which of the groups parameters were used (e.g.: Rp H or Rp W ; 



■FH -FW j. 

% or i , etc. 



We per f orm t wo separate fits using parameters obtained by 



Hebrard et al. 



Winn et al 



(|2010| ) and 



(120091 ) as fixed constraints. These parameters, which we hold constant while 



fitting, are listed in Table [2j Our two sets of results (corresponding to the two different sets 
of fixed parameters) for the planetary ra dius, inclination, and mid- transit time are lis ted 



Hebrard et al. 



(2010) and 



Winn et al. 



(120091 ) for 



in Table |3] along with the results of both 
comparison. We also report the Xr of each of our fits. However, the difference between the 



two yz is minima 



the 



Hebrard et al. 



This leads us to conclude that while our results determined by using 
(120101 ) constant parameters yields a smaller Xr-> t ne accuracy of our two 



fits is statistically comparable. This agreement is demonstrated in Figure [2] as we plot 
our two best f i ts aga inst each other. The solid blue line represents our best fit employing 



Hebrard et al.1 (120101 ) constraints while the red dashed line is our best fit using IWinn et al 
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( 120091 ) parameters. 



We assume that the R p /R* ratio should be the same for both our fits as the transit depth 
of our data does not change. We find that our (Rp/ R*) FW and (R p /R*) FH measurements 
are separated by 1.2cr. Regardless of the difference in our stellar r adius constrain t , we 



still achieve statistically similar results. Our i? p /i?*results straddle 



Hebrard et al. 



fl2010h 



result. We a 



Winn et al 



so fin d that our (R p / R*) is the smallest ratio of the set of results whereas 



(120091 ) determined the largest R p /R* for HD80606b with a difference of 3a. 



The orbital inclination is directly dependent on the stellar radius, making it difficult 
to compare our two i FH and i FW results. We find it interesting that both our i FH and 



FW results are exactly la larger than the inclinations determined by iHebrard et al. 



and 



Winn et al. 



feoioh 



(120091 ) respectively. This is likely due from the blackouts occurring in the 



ingress and egress of our data, which cause a difficulty in fitting the exact curvature of the 
transit. We do not report an independent estimation of the impact parameter, as we do not 
fit for this variable. Instead, we calculate b using Equation [T] to allow for comparison and 
find that our b FH is slightly larger than our b FW but less than la difference. 

Our mid-transit times in both of these complete fits differ from the timing only fits 
listed in Table [T] because they are the best fit times based on the other fitted parameters. 
As we fit for multiple variables, the resulting errors on these times are higher than the 
ones listed above. We still find a mid-transit time earlier than predicted; however, the 
spread between our arid T™ d results is a 3 minute difference, even though they have 

a less than la discrepancy. Interestingly, our T^J^ is the earliest of all our mid-transit 
times, including those re ported in Se c tion 14.11 leading us to believe the mid-transit time 



dependence on the other 



Winn et al. 



(120091 ) parameters is not 



variation. Both of our and T™ d results are within la of 
transit time and are at least 1.2a away from 



Winn et al. 



causing the transit t iming 



Hebrard et al. 



(120091 ) 



(boiOf ) mid 
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We verify our mid-transit time by fitting public spectroscopic data from iHebrard et al. 



(120 10[ ) and our photometric data simultaneously. While the quality of our transit leads to 
large uncertainties in most of our parameters, we determine a mid-transit with less than 
la difference to our T^ Y . We further confir m the period, eccentricity and longitude of 



periapsis determined by 



Hebrard et a 



parameters established by lWinn et al 



(2010). Since these do not deviate from the orbital 



(120091 ). we do not consider it necessary to delve in to 
this finding. We also fit for the impact parameter and R p /R*. However, both values as well 
as the uncertainties associed with them are extremely skewed. We therefore conclude that 
fitting our transit for more than four parameters leads to large inaccuracies in our mode l 



and o nly r eport those paramete rs fitted with the constant constraints from 



Winn et al. 



(120091 ) and 



Hebrard et al. 



(120101 ). 





Constant Parameters 




R* 


M, 


P 


a 


e 


u 


Hebrard et al. 2010 Parameters Fif 


1.007 Rs m 


l.OXMsun 


111.4367 days 


0.455 AU 


0.9330 


300.77° 


Winn et al. 2009 Parameters Flv 


0.968R S un 


XS&Msun 


111.4374 days 


0.4614 AU 


0.93286 


300.83° 



Table 2 : The two groups of p arameters that were held constant while fitting. Fitting runs us- 



ing the 



Hebrard et a 
i 



with the IWinn et al 



(2010) parameters are denoted with FH while best-fit results calculated 
(j2009l ) parameters are designated with FW . 



5. Conclusions 



Employing the optical telescope MOST we observe the January 20 10 transit of 



HD80606b. Using the fitting routine detailed in 



Barnes &: Fortneyl (120031 ) and modifying 



it for the high eccentricity of this planet's orbit, we addre ss two separate goals. We first 



investigate a mid-transit timing discrepancy discovered by 



Hebrard et al 



that our mid-transit time of 2455210.6449 ±0.0034 HJD is within la of 



(2010) and find 



Hebrard et al. 
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1.010 - 




6 8 10 12 14 

Time HJD +2455200 



Fig. 1. — Here we show the MOST photometry full data set, including errors. The photo- 
metric data are binned to one point every 15 minutes . 





Fitted Parameters 




Rp 


Rp/R, 


i 


Mid- Transit Time 


Xr 


ThisWork FH 


0.987lQ ( ggj-R./ U p 


1007+ 00013 

u - iuu ' -0.0062 


89.283° ±0.024 


2455210.6461 ± 0.00781i^JD 


1.01389 


Hebrard et al. 2010 


0.981 ± 0.023R, /up 


0.1001 ±0.006 


89.269° ±0.018 


2455210.6420 ± 0.0010H JD 




ThisWork FW 


0.905 ± 0.032R Jup 


0.0961 ± 0.0034 


89.346° ± 0.022 


2455210.6439 ± 0.00773HJD 


1.10136 


Winn et al. 2009 


0.974 ± 0.030R,/ up 


0.1033 ±0.0011 


89.324° ± 0.029 


2455210.6590 ± 0.0050HJD * 





*Predicted result 



Table 3: This Work F/f represents best fit values using Hebrard et al. 2010 constraints. This 
Work Flv shows best fit values using Winn et al. 2009 parameters. The variables determined 
by the two other groups are shown for comparison. Xr 2 is shown for a statistical comparison 
of the two fits. 
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Fig. 2. — Unbinned close-up of the MOST transit data. The two lines represent our two best 
fits. The blue line is our best fits constrained by Hebrard et al. 2010 parameters, whereas 
the red dashed line is the best fit constrained by Winn et al. 2009 parameters. 
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by 



(2010 ) results and 20 minutes earlier (3a difference) than the predicted value determined 



Winn et al. 



(120091 ) 



As 



Wu k Murravl (120031 ) have 



theoretically ruled out t 



planet causing these variations and 



re possibility for a separate 



Barnes fc O'Brien! ( 120021 ) discarding the hypothesis for 



the existence of an exomoon, we think the transit timing variation is caused by systematic 
errors depending on the contrasting methods of ground-based and space-based observations. 

We also present an independent set of parameters for the planetary radius, orbital 



inclination, and mid-transit time of HD80606b derived from the MOS 
use th e orbital and stellar parameters determined by 



Winn et al. 



pho tometry. We 



((20091) and 



Hebrard et al 



(120101 ) and hold them constant during the fitting process. We find that our two sets 



of results are statistically similar to each other and there are no unu sual discrepancies 



Hebrard et al. 



(120 lOl ) parameters, we 



occurring even with different constraints. With 

further determine that all of our results are within lcr of their results. However, we also 
discover that our fitted planetary radius is 2.2<r smaller than 



Winn et al 



(120091 ) result. 



When fitting for our own independent mid-transit time, we find that our mid-transit 



times differ between each other by 3 minutes but remain 



they are 18 minutes earlier than predicted by 



Winn et al. 



e ss tha n la apart. Furthermore 



(120091 ). We also conclude that 



th e mid-transi t timin g issue is not due to the differing in the orbital or stellar parameters 



of 



Winn et al. 



( 120091 ) compared to iHebrard et al 



(120 101 ). 



As we are currently relying on the observations of only two transits, we lack the 
data necessary to draw any sound conclusions regarding the discrepancies in our results. 
Future observations, both with ground-based and space-based telescopes, are a necessity 
in order to better address the transit timing issue. These observations will also further 
constrain the orbital, stellar, and planetary parameters of HD80606b and aid in bettering 
our understanding of the formation and evolution of this unique system. 
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This Work"' 
Hebrard et al. 2010 
Shporer et al. 2010 





— Hebrard et al. 2010 Parameters 

— Winn et al. 2009 Parameters 







This Work'" - ' • ' 

Winn et al. 2009* - I » - 

! ! 

0.640 0.645 0.650 0.655 0.660 0.665 

Time HJD +2455210 



Fig. 3. — Five different mid-transit timing results from different groups with errors included. 
The top three results in blue were calculated using Hebrard et al. 2010 parameters. The 
bottom two results in red were calculated using Winn et al. 2009 parameters. 
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